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We report on a theoretical analysis, consisting of both numerical and analytic work, of the 
stability of synchronization of a ladder array of Josephson junctions under the influence of current 
induced magnetic fields. Surprisingly, we find that as the ratio of the mutual to self inductance of 
the cells of the array is increased a region of unstable behavior occurs followed by reentrant stable 
synchronization. Analytic work tells us that in order to understand fully the cause of the observed 
instabilities the behavior of the vertical junctions, sometimes ignored in analytic analyses of ladder 
arrays, must be taken into account. 

PACS No.: 74.50.+r, 05.45.Xt, 05.45.-a 

Ladder arrays of Josephson junctions are intriguing systems for a wealth of reasons: the possibility of phase locking, 
or synchronizing, a maximal subset of junctions suggests their use as microwave sources Q; they offer rich dynamical 
behavior, accessible to both theorists and experimentalists, in the field of coupled nonlinear oscillators (with recent 
interest in the prediction and observation of discrete rotobreathers) Q ; their complexity is between that of better 
understood one dimensional serial and parallel arrays and full 2D arrays (e. g. square arrays) and thus they offer a 
nice link between the two geometries; and ladder arrays can, under circumstances that are partially understood, be 
modeled by the discrete sine-Gordon (DSG) equation ||, which is itself a source of research interest among many 
With a desire to understand better the conditions under which stable synchronization can occur, we consider 
ladder arrays biased with uniform dc bias currents 1b greater than the critical currents of the junctions and include 
the effects of current induced magnetic fields (CIMF's) via self and mutual inductances of the cells of the array, (see 
Fig. [I]) Since the array is current-biased above the critical current there will be a nonzero voltage across some subset 
of junctions in the array. These "active" junctions are synchronized, or phase locked, if their voltages, after some 
initial transients, are identical functions of time. Furthermore, previous workers have established that the effects of 
CIMF's may be important in determining the static and dynamics properties of arrays, and so it is natural to 
consider the effects of CIMF's on synchronization as well. 

In this Letter we present numerical and analytical evidence that mutual inductance between cells of a ladder array 
can lead to rich dynamical behavior, including destabilization of synchronization and reentrant synchronization as the 
relative size of the mutual to self inductance of the cells is increased. The instability results for a finite range of values 
of the mutual inductance and occurs when, for the resistively and capacitively shunted junction (RCSJ) model, the 
numerical solutions to the model equations diverge with time. In addition to our numerical solutions we also investigate 
the dynamics analytically. Namely, the two coupled RCSJ equations for the horizontal and vertical junctions can be 
approximately decoupled. The simplified equations allow us to calculate analytically the Floquet exponents, which 
measure the degree of stability of the synchronization. Then, by direct comparison of the analytic results for the 
Floquet exponents, based on the simplified equations, and the numerical results for the Floquet exponents, based on 
the full RCSJ equations, we learn valuable information about the roles of the horizontal and vertical junctions in the 
array. This technique appears to be a powerful way to analyze the relative importance of subsets of junctions that 
experience different local conditions. 

Josephson junctions consist of superconducting islands separated by a thin layer of nonsuperconducting material. 
In the superconductors, the coherent motion of the paired electrons, or Cooper pairs, leads to a wavefunction of the 
form \P = \^>\e l6 , where 9 is the macroscopic phase of the superconductor. The equations describing the dynamics 
of a single junction depend on the gauge-invariant phase difference, or Josephson phase, across the junction, (ft — 
9i — O2 — (27r/<I>o) Ji A • dl, where A, the vector potential due to an external magnetic field, is integrated along a path 
from one side of the junction to the other. $0 = h/2e is the magnetic flux quantum, where h is Planck's constant 
divided by 2ir and e is the electronic charge. We assume in this work that A = 0. 

Consider a ladder array of underdamped junctions with N cells and periodic boundary conditions, as shown in 
Fig. [|. Each junction has a McCumber parameter [3 C = 2irCI cx R 2 /$o, where C(R) is the junction's capacitance 
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(resistance), and I cx is the critical current of a "horizontal" junction. We assume each cell of four junctions is 
described by a self- inductance L and also has a mutual inductance — M (M > 0) with each of its two nearest 
neighbors, where M < L. As an adjustable parameter in the theory, M/L is allowed to range between zero and one. 
Basic physics arguments for a simple ladder with nearest neighbor mutual inductance precludes M/L from exceeding 
0.5 Nevertheless, it is informative and worthwhile to study the behavior of the model equations over the range 
< M/L < 1. Furthermore, such theoretical work could pique the interest of others to try, using modern fabrication 
techniques, to enhance the mutual over the self inductance of the cells and thereby probe the broader range of M/L 
values studied here. Alternatively, one could point to these results for M/L > 0.5 as an interesting portent to the 
behavior that might be observed in arrays with overlapping cells or in three dimensional arrays. 

An analysis based on the RCSJ equations and including the effect of induced magnetic flux leads to a pair of coupled 
equations for the Josephson phases of the horizontal and vertical junctions, 

&0" + 0' + sin + -J-Z Tr ■ XT 1 ■ (Z ■ + 2-0) = 0, (1) 

PL 



/3 c ip" + ip' + asinip + ttX- 1 ■ (Z ■ + 2ip) = 0, (2) 
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where (3l = 2ttLI cx /^q is the dimensionless self-inductance of a cell, and a = I cy /Icx is the critical current anisotropy. 
The prime symbols denote differentiation with respect to dimensionless time, r = t/t c where t c = hbar/2eI cx R. These 
equations are compactly represented in matrix notation, where and ip are A^-dimensional vectors representing the 
the Josephson phases of the horizontal and vertical junctions, respectively. Z_ is an N x N matrix that depends on 
geometry, and X_ is the dimensionless inductance matrix, also N x TV in size. The diagonal terms of X_ represent the 
self-inductance of a given cell, i.e. Xjj = +1, while the mutual inductance of the nearest neighbor cells is represented 
by the terms Xj,j±i = —^l, where fiL = M/L is the dimensionless mutual inductance. All other elements of X_ 
are zero. Eqs. |]and| were solved numerically for 0, ip, 0', and ip' via a fourth-order Runge-Kutta algorithm as a 
function of the parameters N, (3 C , (3l, V>l, and is = Ib/Icx (the dimensionless bias current). The starting configuration 
consisted of random voltages and zero Josephson phases. 

As described elsewhere for a similar system [|J a stability analysis of the solutions to Eqs. [I] and |^ follows by letting 
= 0o + f] and ip = ipo + <5, where 0o aud ipo are solutions to these equations. Eqs. ^ and are linearized with respect 
to 77 and S. The time dependence of the perturbations has the form 77 ~ e xtcT and 8 ~ e AtcT , where A (A) represent the 
Floquct exponents for the horizontal (vertical) junctions. If Re(Ai c ) > or Re(A/j c ) > we expect unstable behavior 
of the array for the given set of circuit parameters. In fact, we are interested in that exponent whose magnitude is 
closest to zero, as that describes the stability of the longest-lived mode of the array. 

Figure §(a) shows the minimum Floquet exponent as a function of fiL for a 5-cell ladder. The three different sets of 
symbols represent three different sets of initial values for the voltages, 0' and ip' ■ The data show that the exponents 
do depend, to some degree, on the values of the starting voltages, although all the runs follow the same trend. As 
fiL is increased from zero towards 0.5, the stability of phase-locking increases, as shown by the negative exponent of 
growing magnitude, while the degree of stability decreases for increasing /xl greater than approximately 0.6. Even 
more interesting is the behavior of the ladder in the range 0.5 < /Jl < 0.6. For these values of the mutual inductance 
the ladder is actually unstable! This is evidenced by very rapidly growing phases and voltages with time as Eqs. |l| 

and H arc numerically integrated. For N — 5, the lower limit of this instability region is /z^ 1 ' = 0.5 independent of 

(2) 

other circuit parameters such as p c and /3l- The upper limit of this region, which we denote by fi L , depends on such 

quantities as the value of the starting voltages as well as on the value of /3l- For example, for a fixed set of starting 

(2) 

voltages, we find that fj, L ' is a decreasing function of (3l- Also, it is interesting to note that this instability region 
does not appear at all if both the phases and the voltages are initialized to zero! (See discussion below for the reason 
for this behavior.) 

Physically, this instability is due to a competition between the self- inductance of a given loop (say loop j), which 
wishes to have a current with a given sense of circulation, and the mutual inductance of the two neighboring loops 
(j±l), which wish to have the current in loop j flow in the opposite sense. We have also looked at ladders with N — 6, 
7, 8, and 9. All show this instability in the vicinity of hl = 0.5. Indeed, we would expect this competition-induced 
instability to be independent of ladder size for the case of nearest-neighbor mutual inductance in that the onset of 
the instability should always occur at fi^ = 0.5. 
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Interestingly, ladders with N — 7, 8, and 9 also exhibit a second instability region that the 5-cell ladder does not 
exhibit. Fig. ||(b) shows the Floquet exponents for N = 7. This second instability region has an onset at a value of 
> (J,j* that is dependent on ladder size. (See cross-hatched region in Fig. §(b).) We now turn to an analytic 
calculation of the Floquet exponents, which helps us understand the source of these instabilities. 

A reasonable starting approximation is to ignore the coupling between the horizontal and vertical junctions but 
otherwise not to ignore the effects of the vertical junctions. That is, we let ip = in eq. [I] and <j> = in eq. |^. An 
analysis like that described in Ref. B applied to eq. [I] leads to a Mathieu equation describing the time dependence 
of the perturbations to the horizontal Josephson phases. In such a case, the corresponding Floquet exponents for the 
horizontal junctions can be calculated analytically. The result is 



Re(A m i c ) = --^±-±-RcJl-4c^ ) (^), (>!) 



2/3 c 2/3 c y \(3 L 

where we can think of u)^ as the effective normal-mode frequencies of the ladder. For example, = 
(4sin 2 (m7r/A0 + 2^ L {cos(2nm/N) - cos(47rm/iV)})/(/4 - fi L - 1), where < m < N - 1. Note that w^ V) is a 
function of fiL. 

Eq. U was used to produce the solid curves in Figs. ||(a) and (b). Note that if, for particular values of u>m \ Pc and 
Pl, the argument of the square root in Eq. ||is positive and larger than one, then at least one of the Floquet exponents 
will be positive, signaling unstable phase-locked solutions. In fact, since f3 c > and (3l > such an instability will 
occur ifw^f < 0! Plots of lu[ NS> versus \il for N = 5 and 7 (not shown here) demonstrate that aj[ 5 ^ > for < /i_L < 1, 
but oj^ is negative for > 0.8. We have checked that Um > for to / 1 and for N — 5 and 7. Thus the cause of 
this instability in the 7-cell ladder for (ii, > /j}^' is the m = 1 normal mode. That is, this instability is a geometrical 
effect, in that it does not occur for N = 5 for example, and it is triggered by an effective normal mode frequency of 
the horizontal junctions becoming negative. 

This analytic work, however, does not point to the horizontal junctions as the cause of the instability near /x^ = 0.5. 
To appreciate this behavior it is crucial to look to the vertical junctions. A procedure similar to that which led to 
eq. H leads to a set of effective Floquet exponents for the vertical junctions 
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where the geometrical factor jrnP is also a function of /i£ and is similar but not identical to uj^\ In this case, the 

vertical junctions will exhibit an exponentially growing Josephson phase if < —{a(3L cosif>o)/2. Now a plot of 7^ 

versus hl (see inset of Fig ^|(a)) shows that the function abruptly becomes negative at /i£ = 0.5 and asymptotically 

approaches zero from the negative side as fxi is increased further. (We have checked that > for m =/= 0. 

Also, we see similar behavior for the 7-cell ladders.) If we assume that cos^o > 0, then the vertical junctions will 

be unstable for 7^ < 0. Based on the behavior of 7^ an instability region will exist for a range of values, 

/i^ 1 < hl < t 1 ^ where fi^ — 0.5 and /i^ 2 ' will depend on a, (3l, and cosf/'o- For example, as (3l increases we expect 
(2) 

that ir L ' will decrease, i.e. approach a value of 0.5. This agrees with the behavior of the numerical results for the 

Floquet exponents. Also, the inequality 7^^ < — (a/?L cos-0o)/2 suggests that the value of fif^ should depend on the 

value of cos • This is relevant to the numerical results in Fig. |^, where we see that the value of fjr£^ does indeed on 
the choice of the starting configuration of phases and voltages. In general, then, it is clear that the instability near 
Hl = 0.5 originates with the vertical junctions and would thus be missed by an analysis that was based solely on the 
horizontal junctions. It is also clear why this instability does not appear numerically when both the Josephson phases 
and the voltages across the junctions are initialized to zero. In such a scenario, although the horizontal junctions 
may be active, the only possible solution for the vertical junctions is to keep zero voltages and Josephson phases for 
all times. Since we know this instability region is triggered by the vertical junctions, the vertical junctions have no 
chance to "go unstable" and thus the instability never appears. 

We conclude that mutual inductance between cells of an underdamped ladder array has the effect of destabilizing 
synchronization for ranges of values of hl, the ratio of the mutual to self inductance. These specific ranges of [1l that 
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lead to unstable behavior are geometry dependent. An analytic calculation of the Floquet exponents based on the 
horizontal junctions agrees with the numerical exponents, based on the full RCSJ equations, for those values of /ir, 
for which stable phase locking occur. To understand the cause of all the observed instabilities, however, it is crucial 
in the analytic work to consider the behavior of the vertical junctions. 

Although some values of the mutual inductance used in these simulations can not be obtained in simple ladder 
arrays, this work suggests that experimentalists may wish to attempt fabrication of arrays that enhances the mutual 
over the self-inductance, perhaps making it possible to look for the rich dynamical behavior predicted here. Certainly 
researchers working on the problem of coherent emission from Josephson junction arrays should be aware this potential 
for unstable behavior exists. 
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useful conversations. This research was funded by the Howard Hughes Medical Institute Undergraduate Biological 
Sciences Education Program grant #71196-529503 to Ohio Wesleyan University. 
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FIG. 1. Ladder array of Josephson junctions with periodic boundary conditions. The horizontal junctions, along the rungs 
of the ladder, are parallel to the x axis, while the vertical junctions are parallel to the y axis. A dc bias current, Ib, is injected 
at each node on one side and extracted from the opposite side. The Josephson phase for the horizontal (vertical) junction in 
the jth plaquette is <f>j {tpj)- 



FIG. 2. Minimum Floquet exponent for periodic ladders versus the dimensionless, nearest-neighbor mutual inductance. ( 
The results correspond to is = 10, /3 C = 10, and /3l = 100.) The three different sets of symbols represent three different 
starting configurations of voltages across the junctions, while the Josephson phase differences were always initialized to zero. 
The solid line represents an analytic result (Eq. ||) based on the horizontal junctions, (a) N = 5. The analytic result predicts 
stable phase-locked solutions for < /xl < 1. The numerical results exhibit an instability, however, for jit^ < /il < where 
fi^/ — 0.5 and fj,^ is dependent on the starting configuration of phases and voltages, as well as on the value of Pl- This 
instability originates with the vertical junctions. INSERT: geometric quantity 7q 5 ' versus ^l- (see eq. ^ (b) N = 7. In this 
case, the geometry of the ladder leads to an instability for fiL > 0.8 that originates with the horizontal junctions. This "large 
/tti" instability is marked in the figure as a cross-hatched region. The instability due to the vertical junctions near fiL = 0.5 
still exists but is narrower than for the 5-cell ladder. 
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This figure "fig2a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0008278vl 



This figure "fig2b.jpg" is available in "jpg" format from: 
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